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ABSTRACT 

Values  of  four  parameterization  constants  for  cne  ouiK, 
secpnd-order  closure  turbulence  model  of  Garwooo  (1976)  are 
determined  using  ocean  data  from  the  Nortn  Paciric  experi¬ 
ment  (NORPAX)  Trans-Pacific  ship  of  opportunicy  program 
(TRANSPAC)  ana  atmospheric  forcing  from  the  Fleec  Numerical 
Oceanography  Center  (PNOC)  mooei.  The  annual  cycle  or  sea 
surface  temperature  is  simulateo  for  eacn  of  tne  years 
1976-1978.  The  difference  between  the  simulated  ano  cne  an¬ 
alyzed  monthly  change  of  sea  surface  temperature  is  mini¬ 
mized  to  determine  the  optimal  values  of  the  model  con¬ 
stants.  >The  model  is  relatively  insensicive  to  cnanges  in 
the  constants  for  this  set  of  data.  However  simulations  us¬ 
ing  the  values  determined  in  this  stuoy  yieia  results  wnicn 
are  quite^ifferent  from  those  obtained  using  vaiues  for  cne 
modelpdnstants  determined  from  earlier  studies  witn  difrer- 
enb^nitialization  and  forcing  data. 

^  Three  statistical  techniques  are  used  to  determine  the 
sensitivity  of  the  model  to  cnanges  of  the  model  conscanrs. 
Two  techniques r  bootstrapping  ana  pooled  permutation  proce- 
durOf  yield  similar  conclusions.  The  thirOf  a  standard  stu- 


1.  INTRODUCTION 


This  paper  aadresses  chree  interrelatea  copies  wnicn 
are  fundamental  to  simulating  the  chermoaynamics  or  cae  up¬ 
per  ocean  on  diurnal  to  seasonal  time  scares.  Txie  rirsc 
topic  is  the  determination  the  mooel  constants  rcr  cne  sec¬ 
ond-order  closure  bulk  mooel  of  turbulent  oceanic  oouiiaary 
layers  developed  by  Garwood  (1976,  1977).  Tne  constants  are 
determined  for  forcing  oeriveo  from  the  Fleet  Numerical  Oce¬ 
anographic  Center  (FNOC)  atmospneric  primitive  equation  moo- 
el  and  ocean  thermal  structure  information  oocaineo  rrom  cne 
North  Pacific  Experiment  (KORPAX)  Trans-Pacific  ship  of  op¬ 
portunity  XBT  (TRANSPAC)  analysis.  Seconoiy,  enree  statis¬ 
tical  tests  are  used  to  examine  the  variations  in  mooel 
skill  as  a  function  of  the  parameterization  constants  ano 
the  solar  radiation  absorption  model.  Finally,  tne  model 
for  the  absorption  of  solar  radiation  in  tne  upper  ocean, 
which  is  used  in  the  Garwood  model,  is  examined  ano  compared 
with  several  other  models. 

Turbulence  closure  is  necessitated  oy  tne  process  of 
Reynold's  averaging  which  generates  correlations  (moments) 
of  the  fluctuating  variaoles,  resulting  in  more  unxnowns 
than  equations.  The  proolem  can  not  oe  resoiveu  at  a  funda¬ 
mental  level  with  additional  equations  since  each  new  equa¬ 
tion  for  a  correlation  introduces  additional  nigner  oroer 
moments.  The  system  of  equations  must  oe  closed  oy  speciry- 
ing  the  highest  order  of  the  terms  to  oe  retained  and  param¬ 
eterizing  these  correlations  in  terms  of  lower  order  moments 
and  mean  values.  One  may  then  cnaracterize  the  various  tur¬ 
bulence  models  according  to  the  oroer  at  wnich  tnis  closure 
is  made  and  the  particular  parameterizations  used.  Tne  Gar¬ 
wood  (1976,1977)  model  discussed  in  the  paper  employs  a  sec¬ 
ond-order  closure  scheme.  Thus,  the  seconu-oruer  moments 
(triple  correlations)  are  parameterized  in  terms  of  first- 
order  moments  (autocorrelations  and  cross  correlations)  ano 
zero-order  moments  (mean  values) . 

The  parameterization  of  the  nignest  order  moments  gen¬ 
erates  constants  of  proportionality  between  the  hignest  or- 
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der  moments  and  functions  of  the  lower  oroec  moments.  Tnese 
constants  must  be  determined  experimentally  before  cne  mouej. 
can  be  used  to  simulate  accurately  the  fluid  flows  of  inter¬ 
est.  The  goal  in  turbulence  modeling  is  to  formulate  tne 
parameterizations  such  that  they  are  "universal"  or  applica¬ 
ble  to  a  wide  range  of  geophysical  flows.  Unfortunately r  no 
parameterization  can  be  completely  universal.  Parameteriza¬ 
tions  derived  for  bouncary  layer  flows,  for  example,  can  not 
be  expected  to  work  well  for  free  mixing  layers,  ]ets  or 
wakes,  since  the  physics  ana  the  bounoary  conoitions  are 
quite  different.  Further,  for  a  given  parameterization,  tne 
optimal  values  of  the  constants  may  also  be  expecteo  to  oir- 
fer  for  different  oceanic  and  atmospheric  regimes  since  ali 
possible  physical  processes  are  not  includeo  in  the  mooel. 
In  addition,  the  space  or  time  resolution  of  the  oata  useo 
may  not  be  adequate.  The  constants  for  this  mooel  nave  pre¬ 
viously  been  estimated  form  laboratory  data  (Garwood,  1976) 
and  ocean  weather  ship  observations  at  Ocean  Station  Papa 
14^)  (Garwood  and  Adamec,  1982)  .  Although  these  ear¬ 
lier  values  were  not  optimal  determinations,  they  will  oe 
compared  with  the  reults  obtained  here  using  FNOC  anu 
TRANSPAC  data. 

The  absorption  of  solar  radiation  in  the  upper  ocean 
significantly  influences  the  diurnal  cycle  of  the  mixea 
layer,  the  summer  minimum  mixed  layer  depth  ana  associatea 
temperature  maximum  and  the  spring  transition.  The  physics 
of  the  absorption  of  solar  radiation  is  well  known;  however, 
the  distribution  of  yellow  matter  and  particulate  matcer, 
which  significantly  affects  the  process,  is  not.  Further, 
an  elaborate  bandwidth  dependent  model  of  solar  raaiacion 
absorption  is  very  computationally  expensive.  Thus  ic  is 
necessary  to  develop  an  accurate  but  simplified  parameteri¬ 
zation  of  the  absorption  process.  The  formulation  used  in 
this  model  will  be  compared  with  other  methods,  incluaihg  an 
empirical  model  (Zaneveld  and  Spinrad,  1980)  which  best  rep¬ 
resents  the  limited  data  sets  that  are  available. 


The  choice  ana  application  of  statistics  wnicn  are  ap¬ 
propriate  to  this  type  of  stuoy  is  a  non-triviai  aspect  or 
the  work.  Three  relatively  simple  statistics  nave  oeen  cnu- 
sen:  one  Gaussian  statistic  ana  two  non-Gaussian  statistics 
which  are  capable  of  generating  their  own  prooaoiiity  uis- 
tribution  functions  (paf's).  It  is  nopeo  that  tnis  latter 
approach  will  reduce  the  problems  causea  by  using  sample 
populations  which  do  not  necessarily  have  Gaussian  paf's  ana 
whose  members  are  not  indepenaent. 

The  Garwood  moael  is  reviewed  in  section  2.  Tne  ocean 
thermal  structure,  surface  forcing  ana  tne  aecails  of  the 
numerical  simulations  for  this  study  are  aiscussea  in  sec¬ 
tion  3.  The  results  of  the  numerical  experiments  are  aa- 
dressed  in  section  4  and  a  summary  is  proviaea  in  section  5. 


2.  DESCRIPTION  OF  MODEL 

Since  the  moael  has  been  describea  in  aetail  oy  Garwooa 
(1976  and  1977)  ,  we  will  only  briefly  review  the  equations 
to  indicate  the  model  constants  which  riust  oe  aeterminea. 

Ga'rwood  integrated  the  equations  of  motion  ana  the 
buoyancy  equation  through  the  mixea  layer  to  yield 

^U> 

hw - 

dt 

^<V> 

hr - 

dt 


-  f  h<V>  =  -uw'(O)  +  uV(-h) 


+  f  h<U>  =  -vw(0)  +  vw(-h) 


(ii.l) 


(2.2) 


\-ti)  -  I  Q - ij,  j  Q  az' laz  . 


-bw(0)  +  bw 


The  vertical  mean  value  of  any  variaole,  X,  is  uerinea  as 


(2.3) 


<X> 


X(z)  dz 


where  h  is  the  mixed  layer  depth  aha  z  is  zero  at  the  sea 
surface  and  positive  upward.  The  boundary  conditions  at  z  = 
0  and  z  »  -h  are 


-uw^(  0) 


uV(-h)  =  -^A_  Au  ““ 
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-  vV(  0 )  =  / 
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vw*{-h) 


-_A_  Av  —  / 


bw*{0) 


ocqQy  ;ig(E-P)S 


uc 

un 


Dw'(-n)  =  -_A_  A^  —  • 

ac 


and  ^  are  the  x  and  y  components  of  the  surface  stress. 


is  the  net  surface  heat  flux  (total  heat  flux, 


less 


solar  radiation^  Q^) .  S  is  the  surface  salinity  ano  (E-P) 


is  the  evaporation  less  precipitation  race.  cne  Heav¬ 

iside  step  function,  is  oefinea  as 


ah 


_y\_  =  0  for  —  <  0 
dt 


on 


and 


_A_  =  1  for  —  >  0  . 

at 


The  changes  in  mean  velocity  components  ana  ouoyancy  tnrougn 
the  entrainment  zone  are 

Au  =  U(-h)  -  U(-h  -  S)  , 


Av  =  V(-h)  -  V(-h  -  , 

Ab  =  B(-h)  -  B(-h  -  5)  . 

Ao.  Av  and  Ab  are  the  changes  in  the  mean  properties 
across  the  entrainment  zone  and  dh/at  is  the  rate  of  mixea 
layer  deepening. 

The  turbulence  equations  integratea  througn  che  mixea 
layer  are _ 

(AU)^  +  (AV)-*'ldh 


d(h<u^  +  v^> 


dt 


•j 


dh 

-  —  -  m,<e>^ 

Jot  * 


^Ji(<e>  - 


3<v 


-  +  — fhj<'e>  , 


(2. 


d(h<w^>)  h  bV(0)  -  dh  o<>g 

dt  2  2  at  ^ 


+  m. 


—  —  m,  r _ V  ra-1  — 

,<e>3i(<e>  -  3<w^>) - <eA  +  -*fn  <e>  . 

‘  3  . 


(2. 


m. 


J 


Integrating  the  total  turbulent  kinetic  energy  (TKE)  tnrougn 
the  entrainment  zone  provides  the  final  equation  neeaea  co 
close  the  system ^of  equations, 
dh  2^  (wX/k  <e> 

dt  hAs  +  2<e> 

is  the  vertically  averaged  TKE 


(2 


The  left  hand  side  of  (2.4)  is  the  local  rate  of  cnan^e  of 
the  vertically  averaged  horizontal  TKE  in  the  mixea  layer. 

The  time  scale  for  the  terms  on  the  left  nano  siae  of 
(2.3)  and  (2.4)  is  approximately  100  seconds  for  cne  oceanic 
planetary  boundary  layer  unoer  normal  circumstances  (ae 
Szoeke  and  Rhines»  1976).  Since  we  are  interescea  in  montn- 
ly  and  longer  time  scales  in  this  stuoy,  tnese  terms  will  oe 
neglected.  The  first  term  on  the  rignt  nano  siae  of  f'^.3) 
is  a  parameterization  of  the  surface  winu  stress  actj  on 

the  shear  of  the  mean  current,  the  proauction  cue  to  o  k- 

ing  surface  waves  and  the  transport  of  horizontal  TKE. 
is  the  ocean  friction  velocity,  given  ^  ”  i^l/  ^  i 

used  as  a  scale  for  the  rms  velocity  fluctuations.  Tne  sec¬ 
ond  term  on  the  right  hand  side  of  (2.3)  is  the  momentum 
flux  associated  with  the  entrainment  of  nonturouient  fluia 
at  the  base  of  the  mixed  layer.  The  momentum  flux  uue  to 
internal  waves  is  neglected.  The  third  term  is  a  parameter¬ 
ization  of  the  pressure  rate-of-strain  terms  wnich  act  to 
reduce  the  anisotropy  of  the  turbulent  flow  (Rotta,1951; 


Lumley  and  Khajeh-Nouri,1974) .  The  last  term  is  a  parame¬ 
terization  of  the  dissipation  rate.  The  Dissipation  rate 
can  be  estimated  from  the  rate  at  which  the  large  scale  ea- 
dies  supply  energy  to  the  smaller  scales.  As  aiscussea  in 
Tennekes  and  Lumley  (1972) ,  the  dissipation  rate  is 
du;  ^u; 

^  *9  9* 

r 


E 


C  «U;  du: 

i  r  O  y  dz  «  — 
J  dxj  dx^  T^ 


where  T^^  is  the  dissipation  time  scale.  There  are  two  time 
scales  for  the  oceanic  mixed  layer:  the  time  scale  imposea 
by  planetary  rotation,  ana  the  time  scale  rormea 
from  the  turbulent  velocity  and  the  length  scale  of  the 


large  scale  turbulent  flow,  T^*^  h/u. 


The  first  scale  is 


"r  **'  ♦* 

the  characteristic  time  scale  of  vortex  stretching  aue  to 
the  planetary  rotation.  The  secona  is  the  time  requirea  for 
a  typical  large  scale  eody  to  overturn.  A  feature  of  Gar¬ 
wood's  model  is  that  both  of  these  time  scales  are  incorpo- 


rated  into  the  parameterization  of  the  oissipacion  ciu.e 
scale  by  specifying  that  is  scalea  wich  the  smaller  of 
or  T^.  To  incorporate  both  time  scales  the  oissipacion  time 

scale  is  defined  as 

•I  -I  -I 

Te  =  Tr  +  Tt  . 

In  (2.5)  ,  the  first  two  terms  on  the  rignt  nano  siae 
are  the  surface  and  entrainment  buoyancy  fluxes.  Tne  cniro 
term  is  the  absorption  of  solar  raoiation  wicn  oepcn.  Tne 
remaining  terms  are  identical  to  those  in  (2.4).  (2.1) 

through  (2.6)  form  a  complete  set  of  equations  for  the  ae- 
pendent  variables  <U>r  <V>,  <B>,  <u*»  +  v>*>,  <v/K>  ana  n. 

The  depth-depenaent  absorption  of  solar  raoiation  i.s 
particularly  important  in  shallowr  summer  mixea  layers.  In 
tnese  cases  the  amount  of  solar  raoiation  aosoroeo  in  tne 
mixed  layer  significantly  affects  the  mixed  layer  tempera¬ 
ture,  while  the  absorption  below  the  mixeo  layer  arreccs  tne 
entrainment  rate  and,  therefore,  the  mixeo  layer  oeptn.  Tne 
flux  of  solar  radiation  in  the  ocean  can  b'  written  as 
Qjz)  «  Qj<*)^ 

where  A  is  the  wavelength  ana  v;(A)  is  the  optical  water  type 
(Jerlov,  1976).  Various  parameterizations  (see  Table  1) 
have  been  proposed  for  the  adsorption  of  solar  raoiation. 


TABLE  1 

Parameterizations  of  solar  radiation  in  tne  ocean 


Parameterizations 


Q^z)  = 

Qj(z)  =  (^•)(r,  e' 
qJz)  =  r;  e^*** 

Q^z)  -  (^j<|e^*(l  -  a  tan’*  (ox)) 
(^z)  »  Q<iKl-r) 


+  i:,e*'»*) 


Izl  <  1 
Izl  >1 


Source 
Denman,  1973 

Paulson  and  Simpson,  1977 
Konoo,  1979 

Zanevelo  ano  Spinrao,  19ti0 
Garwooo,  1976 


Q.«ii)is  the  solar  raoiation  at  the  sea  surface,  tne  s  are  tne 
extinction  coefficients  and  the  r's  are  constants  (0  <  r  <  1) . 


The  first  three  of  these  parameterizations  assume  cnat  a 
significant  fraction  of  the  solar  radiation  absorption  can 
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be  modeled  by  using  one,  two  or  many  wavelength  oanos,  re¬ 
spectively.  The  fourth  parameterization  was  ootainea  eui^ir- 
ically  using  an  extensive  data  set  of  absorptance  anu  crans- 
mittance  measurements  (Spinrad  et.  ai. ,  1979).  Tne  con¬ 
stants  in  the  fourth  formula  in  Taole  1  oepeno  on  cne  water 
type.  For  this  study  the  appropriate  water  type  is  Type  11 
(Zaneveldr  personal  communication) t  and  the  constants  are  M 
3  0.067  m/Sr  a  »  0.4158,  and  b  »  3.9865  m/s.  Simpson  anc 
Dickey  (1981)  comparec  the  first  four  parameterizatxons  in 
Table  1  and  founa  that  using  a  single  exponencial  mouei  is 
acceptable  only  at  very  high  winu  speeas  (U,^  >  20  m/s) . 
They  also  founa  that  the  empirical  formulation  of  Zaneveiu 
and  Spinrad  (1980)  and  the  double  exponential  mooel  or  Simp¬ 
son  and  Paulson  (1977)  produce  nearly  identical  temperacure 
profiles,  and  that  the  complexity  of  the  spectral  aecomposi- 
tion  form  is  not  necessary  in  mixed  layer  modeling. 

The  fifth  parameterization,  used  in  the  Garwooa  moael, 
is  a  double  exponential  model  for  which  one  of  the  extinc¬ 
tion  depths  approaches  zero.  A  fixed  fraction  (1-r)  of  Q 
is  absorbed  in  the  first  meter  and  the  remaining  flux  is  ab¬ 
sorbed  exponentially  for  |z|  >  Im.  Values  of  r  ana  o  must 
be  specified.  This  parameterization  will  be  comparea  with 
Zaneveld's  empirical  formulation. 

In  (2.4)  to  (2.6)  ,  m,  through  are  the  constants 
which  must  be  determined  from  laboratory  and  geophysical 
data.  Of  these  only  four  are  inaepenaent,  ana  only  three  (ir^ 
and  p^  -  m^/m,  ) 

are  critical  to  moaeling  the  annual  cycle  of  the  upper 
ocean  mixed  layer.  In  this  study,  we  use  hear  ana  momencum 
fluxes  derived  from  the  FNOC  atmospheric  fielas  ana  cne 
TRANSPAC  ocean  data  over  a  range  of  latitudes  ana  longitudes 
to  estimate  the  optimal  values  for  m^,  p^»  ,  r  ana  X, 

3.  DESCRIPTION  OF  NUMERICAL  EXPERIMENTS 

Choosing  an  annual  time  periou  for  this  stuuy  requires 
that  the  moael  be  able  to  simulate  four  significant  pnenome- 
na  each  of  which  requires  different  responses  from  tne  moo- 
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el.  The  fall  deepening  is  pe  'ticularly  sensitive  co  ruj  (eg. 
2,4)/  the  parameterization  of  surface  wino  mixing,  Tne  win¬ 
ter  minimum  mixed  layer  temperature  is  sensitive  zo  (egs. 
2, 4, 2, 5)/  tne  dissipation  paramecer ization.  The  spring 
transition  depends  on  the  interaction  of  m  ano  tne  surface 
heat  fluxes.  Finally,  the  summer  maximum  mixeo  layer  temp¬ 
erature  depends  on  r  and  ^  in  the  parameterization  of  solar 
radiation  absorption, 

A  grid  covering  a  large  area  of  the  North  Pacific  rroiu 
34*H  to  46*N  in  4*  increments  ano  from  170*W  co  140* W  in  10* 
increments  was  chosen  for  this  scuoy.  This  area  is  large 
enough  to  encompass  various  oceanic  tnermal  scrucrure  ana 
atmospheric  forcing  regimes  but  ooes  not  incluoe  areas  wnicn 
are  strongly  influenced  by  eastern  or  western  bounoary  cur¬ 
rents  where  horizontal  processes  (e,g,  aovection  ana  aiver- 
gence)  are  likely  to  be  important.  The  north-south  extent 
of  the  grid  incorporates  only  the  area  whicn  nas  consiscenc 
data  coverage  throughout  the  entire  year.  Regions  to  cne 
north  (south)  of  these  latitudes  have  fewer  data  points  cur¬ 
ing  the  winter  (summer)  due  to  the  latitudinal  shifc  of  snip 
tracks  with  the  seasons. 

The  specific  quantity  chosen  to  compare  the  mouel  re¬ 
sults  wich  the  data  is  the  monthly  change  in  sea  surface 
temperature  (SST)  averaged  over  an  annual  time  perioa  ano 
over  the  grid.  The  monthly  change  in  SST  is  more  difficult 
to  simulate  than  the  SST  since  the  large  scale  nortn-souch 
SST  gradient  in  the  north  Pacific  contains  more  of  tne  vari¬ 
ance  than  that  associated  with  the  monch-to-raontn  cnange  in 
the  SST  field,  va<-cnticy  co  ot  uiinimizeo  can  oe  'wriccen 


where  M  is  the  number  of  grid  points,  N  is  the  numoer  of 
months  used  in  the  average,  ano  At^^  ana  At.  are  tne  moncnly 
changes  in  the  simulateo  and  analyzeo  SST  respectively. 
This  quantity  was  choosen  rather  than  an  rms  airference 
since  it  allows  the  study  of  some  aadicional  scacisicicai 
techniques  which  will  oe  reportea  elsewnere. 
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The  ocean  thermal  structure  fields  wnicn  are  useu  in 
the  initialization  of  the  model  and  in  tne  comi^utacxon  of 
the  model  error  were  obtained  from  the  NORPAX  TRANSPAC  anal¬ 
ysis  (White  ana  Bernstein,  1979) .  Tne  optimal  incerpoiation 
analysis  uses  the  temperature  ooservations  aurin«j  a  moncn  to 
estimate  the  temperatures  at  each  grid  point.  Horizontal 
analyses  are  prepared,  indepenaently ,  at  0 ,  20,  40,  t>0,  80, 
100,  120,  150,  200,  250,  300,  and  400  mecers.  These  rieids 
are  interpolateo  to  5m  depth  intervals  using  a  conuDinacion 
of  linear  and  exponential  interpolation  (Elsoerry,  ec.  ai. , 
1982) .  The  model  is  initialized  on  January  15th  anu  inte¬ 
grated  through  December  15th  for  1976,  1977  ana  1978.  The 
model  generated  monthly  SST  used  in  (3.1)  is  a  five  aay  av¬ 
erage  centered  on  the  15th  of  each  month. 

The  surface  boundary  conditions  requirea  for  tne  Gar¬ 
wood  mooel  are  the  total  heat  flux  (Q^) ,  the  solar  radiation 
(Q^)  and  the  wind  speed  (W) .  These  variables  are  oeriveu 
from  fields  generated  by  the  FMOC  atmospheric  preaiction  ana 
analysis  system.  Instantaneous  values  of  ana  are  ar¬ 
chived  at  12-h  intervals  and  W  at  6-h  intervals  on  a  63  x  63 
polar  stereographic  grid  of  the  northern  hemisphere.  Bili¬ 
near  interpolation  is  used  to  map  these  fields  onto  a  4*  la¬ 
titude  by  10*  longitude  grid  chosen  for  this  stuuy.  Tne 
values  are  interpolated  to  1-h  intervals  (Gallacner,  1979) 
and  is  corrected  to  remove  a  systematic  long  cerm  bias 
(Elsberry  et.  al. ,  1982). 

The  procedure  used  to  find  one  oest  fic  of  one  mooel 
simulations  to  the  analyzea  fields  was  to  first  set  r  =  0.5 
and  it  -  0.1  ra"*  and  minimize  (3.1)  in  (ir^,  p^)  space  for 
each  year.  Next,  r  was  held  conscant  ana  m^ ,  anu  were 
varied  to  achieve  the  best  fit  to  the  1977  annual  cycle. 
This  value  of  ii  was  comparea  to  the  optimal  founa  inaepen- 
dently  by  comparison  with  Zanevela's  empirical  formula.  In 
the  comparison  with  Zaneveld's  formula  botn  r  anu  if  are  var¬ 
ied. 


TTx."  »- 


4.  DISCUSSION  OF  RESULTS 

As  a  preliminary  test»  che  1977  annual  cycle  was  useu 
to  determine  the  distribution  of  the  model  error  in  (ir^r  ^ 
space.  A  broad  but  physically  reasonable  ranije  or  anu 
was  tested  using  widely  spaced  values  of  .  Tnese  sim¬ 

ulations  indicate  a  local  minimum  of  E  in  the  neignoornooo 
of  m^  s  2.0,  p^  »  4.0.  The  minima  for  all  three  years  were 
located  by  testing  the  region  centereo  on  p^=  4.0,  2,0 

(Figs,  la,  b  and  c) .  For  each  year  the  minimum  lies  on  tne 
line  p^=  4.0,  but  the  value  of  m^  varies  oecween  1.5  ror 

1978  ana  4.6  for  1976. 
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4.1  STAIISTICS  Q£  VARIATIONS  m  XH£  MODEL  CONSTANTS 

The  student's  t-test  is  used  to  determine  wnerner  cne 
variations  in  E  are  statistically  significant.  The  null  ny- 
pothesis,  H^,  is  E(l)  -  E(2)  -  0  where  1  ana  2  represent 

any  two  values  of  the  spatial  and  temporal  average  moael  - 
analysis  differences  defined  by  (3.1),  The  hypothesis  is 
accepted  if 

-t  <  t  <  t 

o(^2,m+n-2  o(/^2,m+n-2 

where  t  is  the  student  variate,  oCis  the  significance  level 
and  m  and  n  are  the  number  of  points  in  the  samples  (m  >  16, 
n  alO  ;  MxN  s  160,  for  these  studies).  The  stuaent  variate 
is  defined  as 

t  =  -  E(2) 

sp  (1/m  +  1/n)  ' 

where  sp,  the  pooled  variance,  is 
(m-1)  si  +  (n-1)  s2 

sp  = -  ■  ■!  -I  -  .  ■■  , 

(m+n-2) 

si  and  s2  are  the  variances  of  E(l)  ana  £(2)  respectively. 

For  the  values  in  1976  (Fig.  la)  ,  the  stuuent's  t-test 
indicates  that  all  non-adjacent  points  are  signir icantly 
different  at  the  90%  confidence  level  (which  will  be  usea 
for  all  cases) .  Furthermore,  tne  oniy  aa3acent  pairs  of  (m^ 
f  Pj)  points  which  do  not  aiffer  significantly  are  (3.0,  4.0 
and  3.5,  4.0),  (3.2,  5.0  and  3.5,  4.0)  ano  (3.2,  5.0  anu 


3.2r  4.0).  In  the  1977  simulations,  the  resuxcs  are 

less  sensitive  to  variations  of  the  parameters.  Tne  pome 
«  1.5  was  significantly  aifferent  from  the  points  uij  =2.2 
and  iQj  *  2.4  with  4.0.  None  of  cne  other  poinets  were 

significantly  aifferent  for  1977.  In  1978,  wxen  p^  =  4.U, 
the  point  m^  =  1.8  differs  significantly  from  =  1.5,  1.6, 
and  3.0  and  m^  *  2.6  is  significantly  aifferent  from  = 
1.4  and  m^  =  1.5.  If  the  constants  are  set  at  p^  =  4.0, 

=  2.2,  the  moae 1-ana lysis  aifferences  are  signiricant,  at 
the  90%  level,  among  the  three  years.  Furtnermore,  tne  min¬ 
imum  in  the  parameter  space  for  each  year  is  significantly 
different  from  the  minimum  in  each  of  the  other  two  years. 

The  student's  t-test  assumes  that  the  samples  are  arawn 
randomly  from  a  normal  population.  Thus  the  aata  points  are 
assumed  to  be  indepenaent  of  each  other.  This  assumption 
can  be  tested  by  pooling  the  A.T' s  from  two  sets  of  moael 
simulations.  A  probability  distribution  function  (puf)  can 
be  generated  by  randomly  drawing  samples  from  tne  new  popu¬ 
lation.  This  pdf  can  then  be  useo  to  determine  tne  signifi¬ 
cance  levels.  Using  this  method  aoes  not  require  that  the 
samples  have  a  Gaussian  pdf  or  that  the  points  be  inaepen- 
dent.  An  example  of  this  procedure,  the  pool  ana  permuta¬ 
tion  procedure  (PPP)  is  describee  by  Preisenaorf er  ana  Moo- 
ley  (1982) .  First,  the  two  data  sets  are  pooled  into  one 
set  where  the  first  m  elements  are  the  first  sample  of  /Vt' s 
and  the  following  n  elements  are  the  secono  sample.  The 
next  step  is  to  permute  the  m  -t-  n  data  points  to  construct  a 
new  pooled  data  set.  The  first  m  elements  become  a  new 
first  sample  ana  the  following  n  elements  become  the  secono 
sample.  The  test  variable  t  may  then  oe  calcuiatea  for 
these  two  new  data  sets.  These  steps  are  now  repeated  for  a 
total  of  100  random  experiments  yielding  100  values  of  t. 
The  t  are  then  ordered  in  increasing  value  witn  t(l)  oeing 
the  smallest  ana  t(lOO)  the  largest.  Then  the  values  of 
t(5)  ana  t(95)  now  yield  the  90%  significance  levels  airect- 
ly,  which  eliminates  the  assumptions  required  to  use  the 
student's  t-test. 
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A  Similar  method  callea  'bootstrapping*  can  oe  usee  zo 
generate  the  pdf  (Preisendorfer ,  personal  communicacioa) . 
In  bootstrapping,  the  selections  are  made  ranoomly  rrom  tne 
entire  data  set,  whereas  only  che  permutations  of  tne  pooieo 
population  are  used  in  the  PPP.  Therefore,  since  repeti¬ 
tions  of  the  data  points  may  enter  into  the  pooieo  sec  using 
bootstrapping,  slightly  broader  reference  oistrioucions  may 
be  generated.  Thus,  the  null  hypochesis  woulo  ceno  co  oe 
accepted  more  often  using  bootstrapping  than  using  tne  PPP. 

The  critical  points  for  accepting  the  null  nypocnesis 
may  be  quite  different  from  che  student's  t-cesc  ir  cne 
original  data  sets  do  not  have  approximately  Gaussian  pars. 
If  this  is  the  case,  fewer  runs  will  be  statiscically  aif- 
ferent  using  the  bootstrapping  methoo  ana  PPP.  Resuics  are 
tabulated  in  Tables  2,  3  ana  4  for  cases  in  wnicn  the  scu- 
dent's  t-test  rejected  the  null  hypothesis  at  the  90%  conri- 
dence  level.  Variations  in  m^  and  within  a  given  year 
(Table  2)  do  not  lead  to  any  significant  differences  except 
using  the  PPP  between  mj  =  2.0  and  m_j  =  3.5  for  p^  =  4.0  in 
1976. 


TABLE  2 

Comparison  of  PPP  and  Bootstrapping 


Year  (m^,  P,)  (nij  rP,) 

1976  (3. 5, 4.0)  (2. 0,4.0)  -11.49 

1976  3. 2, 4.0  i3.2,5.0i  -  3.09 

1976  3. 2, 4.0  2. 2, 4.0  -  8.11 

1977  1, 1.5, 4.0  12.4,4.0  1.86 

1978  1.5, 4.0  3.0,4.0i  -  5.78 


PPP 

t(5)  t(95) 

-11.05  8.66* 

-11.96  10.01 
-12.31  8.02 

-  6.91  8.97 

-  9.46  7.73 


Bootstrapping 
t(5)  cT95) 

-12.58  9.72 

-16.27  10.69 
-12.15  9.62 

-  8.58  7.93 

-  7.65  7.73 


Student  variate  and  90%  confidence  limits  as  determined  oy  tne 
bootstrapping  method  and  by  PPP  for  the  cases  for  wnicn  cne  student's 
t-test  reiected  the  null  hypothesis  at  tne  90%  level  witnin  each 
particular  year.  The  entry  marked  with  an  aster isx  was  tne  oniy 
case  for  which  either  test  determined  that  che  points  were 
significantly  different. 


Holding  the  parameters  constant  (p^  =  4.0,  ra^  =2.2)  ana 
testing  between  years  (Table  3) ,  the  PPP  test  leaas  to  tne 
conclusion  that  the  results  for  1976  and  1977  are  signifi¬ 
cantly  different,  but  tne  bootstrapping  teennique  aid  not. 
Both  procedures  found  1976  ana  1978  to  oe  significantly  air- 
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ferent  for  ra^  »  3.0,  »  4.0  but  not  for  »  2.2, 

(Table  3).  Results  for  1977  ano  1978  are  not  statibcicaliy 
different.  Testing  for  differences  between  mininria  or  eacn 
year  (Table  4) ,  both  tecnniques  founo  that  the  1976  reauica 
are  significantly  different  from  1977  ana  1978,  ouc  1977  ana 
1978  are  not  significantly  different. 


TABLE  3 


Comparison  of  PPP  and  Bootstrapping  Becween  Years 


i 

i 


Years 

t 

PPP 

t(5)  t(95) 

Bootstrapping 
t(5)  t(95) 

vs  1977  1 

[2. 2, 4.0 

1  -  8.01 

-  7.13 

9.38* 

-  8.82 

9.21 

976 

vs  1978  1 

2.2,4.0; 

(  -  4.88 

-  7.84 

7.86 

-  7.75 

7.43 

977 

vs  1978  1 
vs  1978  1 

2. 2, 4.0 
3. 0,4.0 

1  -1.83 

-10.51 

7.58 

-  7.61 

8.81 

976 

1  -13.16 

-12.08 

10.34* 

-  9.53 

8.93* 

Student  variate  and  90%  confidence  limits  as  aeterminea  by  tne 
bootstrapping  method  ana  by  PPP  for  the  cases  tor  wmcn  tne  stuaent's 
t-test  reiected  the  null  nypothesis  at  the  90%  level.  The  entries 
marked, with  asterisks  were  the  only  cases  for  wnicn  either  test 
determined  that  the  points  were  significantly  aifferent. 


TABLE  4 

Comparison  of  PPP  ana  Bootstrapping  Between  Minima 


Years 

PPP 

Bootstrapping 
t(5)  tTIS)^ 

t 

t(5)  t(95) 

1976  vs  1977 

-17.49 

-  6.59 

7.47* 

-12.18  5.86* 

m  Ji  m 

illl 

-  7.42 

-  7,88 

7.44* 

6.00 

-  9.22  6.90* 

-  9.35  6.47 

Student  variate  and  90%  confidence  limits  as  determinea  by  the 
bootstrapping  method  and  by  PPP  for  the  minima  between  eacn  of 
the  three  years.  These  (m-.a  )  minima  are;  1976  (4. 6, 4.0); 

1977  (2. 4, 4.0);  and  197Ml?5 ,4 .0 )  . 

The  entries  marked  with  asterisks  were  cases  for  wmcn  cne  cesc 
determined  that  the  points  were  significantly  aifrerent. 


For  most  cases,  the  bootstrapping  methoa  yieias  smaller 
values  of  t(5%)  and  larger  values  of  t(95%)  than  aoes  cue 
PPP.  In  addition,  both  techniques  yield  smaller  values  of 
t(5%)  and  larger  values  of  t(9S%)  than  the  stuaent's  t-test. 
Again,  this  implies  that  the  bootstrapping  methoa  ana  tne 
PPP  will  find  two  cases  to  be  significantly  aifferent  less 
often  than  will  the  student's  t-test.  This  large  oifrerence 
in  the  critical  values  used  for  the  different  types  of  tests 
indicates  that  the  assumption  of  a  normal  population  is  not 
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valid  and  the  data  points  are  not  indepenoent.  Tnus  cne 
student's  t-test  is  rejecting  when  it  shoulo  oe  accepteo. 

4.2  EXAMINATION  fi£  SOLAR  RADIATION  liQQ£L 

Two  approaches  were  used  to  determine  optimal  vaxues 
for  the  constants r  r  and  in  Garwood's  solar  raaiacion  pa¬ 
rameterization.  In  the  first  methodf  regression  analysis 
was  used  to  determine  the  values  of  r  ana  ^  thac  provioeo 
the  best  agreement  between  Garwood's  ana  Zanevelo's  parame- 
terizations.  Then  r  and  X  were  cnangea  ana  ana  ^  were 
varied  to  minimize  E  for  1977. 

Since  Garwood's  parameterization  is  an  approximacion  to 
a  two  bandwidth  model  of  solar  raaiation  adsorption,  a  non¬ 
linear  regression  was  performed  to  compare  Zaneveld's  param¬ 
eterization  with  that  of  Paulson  ano  Simpson  (see  Taole  1) . 
Zaneveld's  model  of  solar  raaiation  absorption  is  an  empiri¬ 
cal  representation  of  the  extensive  absorptance  ana  crans- 
mittance  data  set  collected  by  Spinrad  et.  al.  (1979);  tnus 
it  is  assumed  that  agreement  with  Zaneveld's  mooel  is  tanta¬ 
mount  to  agreement  with  the  data  set.  The  values  wnicn  gave 
the  best  agreement  in  the  upper  10m  were  (r  =0.4,  - 

O.OSm**)  and  (r  «  0.6  and  ^  =  3.0m'*).  The  assumpcxon  thac 
one  of  the  two  extinction  depths  is  much  smaller  chan  cne 
other  is  verified  in  that  ^  is  approximately  cwo  orders  of 
magnitude  larger  than  ^ .  The  upper  10m  was  cnosen  since 
the  2d)sorption  of  solar  raaiation  above  the  extinction 
depth,  is  most  critical  to  the  mixed  layer  heac  ana  ouoyancy 
budgets.  Linear  regression  analysis  was  also  used  to  com¬ 
pare  the  logarithm  of  Garwood's  parameterizacion  to  che  log¬ 
arithm  of  Zaneveld's  empirical  formulation  for  -10m  <  z  < 
-Im.  That  is,  the  equations 

ln(r)  *  X  z 

and 

Hz  +  ln(l  -  a  tan*(bx)) 

were  compared  to  obtain  the  values  of  r  ana  ^  wnicn  gave  the 
best  agreement.  The  best  agreement  was  obcaineu  for  (r  = 
0.4,  •  O.OSmT*),  Which  again  verifies  thac  Garwooa's  param- 


etecization  approximates  the  double  exponential  mooeir  wnicn 
was  found  to  adequately  describe  the  absorption  of  solar  ra¬ 
diation  with  respect  to  mixed  layer  dynamics  (Simpson  ano 
Dickey,  1981). 


We  are  concerned  with  the  effects  of  the  parameteriza¬ 
tion  of  solar  radiation  on  the  heat  budget  or  the  mixed  lay¬ 
er.  With  this  in  mind,  Zaneveld's  and  Garwood's  parameter i- 
zations  were  compared  by  calculating  the  cnange  in 
temperature  of  an  initially  nomogeneous  water  column  oue  to 
the  absorption  of  solar  radiation.  The  initial  temperature 
of  the  water  column  was  l/c  ano  the  equation 
d  T  1  S  Qiz) 


a  z  * 

was  integrated  for 


12  hours. 


The  surface  solar  raoiacion 


had  a  cosine  dependence  in  time  with  a  peak  at  local  noon. 


The  temperature  profiles  after  12  hours  for  three  cnoices  of 
r  and  o  are  shown  in  Fig.  2a.  The  results  for  (r  »  0.4,  ^  » 
0.08m')  and  (r  =  0.45,  X  =  0.095m'*')  both  overestimate  tne 


temperature  at  Im,  underestimate  the  temperature  at  2m  ano 
slightly  overestimate  the  temperature  from  4  to  10m.  Tne 
set  (r  a  0.5,  ^  a  O.llm"')  estimates  the  temperature  at  Im 
and  2m  well  but  overestimates  the  temperature  from  2  to  10m. 
The  change  in  heat  content  is  best  estimated  oy  (r  =  0.45,  if 
s  0.095m  ) ,  since  the  underestimate  of  the  temperature  in 
the  first  2m  is  approximately  balanced  oy  the  overestimate 
of  the  2  to  10m  temperatures.  Next,  che  value  of  r,  tne 
fraction  of  solar  radiation  whicn  penetrates  1  m,  was  nelo 
constant  and  o  was  varied.  For  r  »  0.5  ano  0.45,  if  »  O.llmT' 
and  o  »  0.095m“',  respectively,  gave  the  oest  agreemenc  oe- 
tween  Zaneveld's  and  Garwood's  parameter izations.  In  Fig. 
2b,  profiles  are  plotted  for  r  =  0.50  ano  if  =  0.10m*', 
0.15m  ,  and  0.20m  .  Fig.  2c  is  similar  to  Fig.  2b,  out  for 
r  •  0.45.  In  both  of  these  figures,  all  three  choices  of  if 
give  good  agreement  near  the  surface  wnere  the  gradient  wicn 
depth  is  very  large.  Compared  with  Zaneveld's  curve,  tne 
curve  for  which  if  ^  0.20m'  turns  too  soon,  wnicn  implies 
that  the  scale  depth  for  absorption  is  too  deep.  Tne  curve 
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for  which  ^  •  0.10m*  has  approximately  equal  area  aoove  ana 
below  Zaneveld's  curve;  which  makes  the  integracea  neacing 
in  the  upper  10m  the  same. 

The  annual  cycle  for  1977  was  simulated  wich  r  =  0.50 
and  ^  •  0.10m*,  0.15m  *  and  0.20m  *  for  a  range  of  values  of 
and  p^.  The  values  of  E  were  higher  for  it  -  O.lSm”*  ana  V  = 
0.20m*'  than  for  if  »  0.10m  *  and  the  minima  were  oisplaceo  to 
higher  m^^  values  for  higher  ^  values.  The  cnange  in  can 
be  explained  since  larger  iC  values  imply  shallower  extinc¬ 
tion  depths.  Therefore,  more  solar  raoiacion  is  aosoroea 
near  the  surface  and  larger  surface  stresses  are  requireo  co 
mix  the  more  stable  water  column.  Using  tne  oootscrapping 
method  or  the  PPP,  there  is  no  significanc  aifference  among 
results  for  V »  0.10m  *,  =  0.15m*  ana  JT  =  0.20m”'  ac  the  90% 

confidence  level,  where  (m^*  2.5,  p^»  4.0). 


4.3  COMPARIfiQM  WIIH  PREVIOUSLY  CALCULATED  VALUES 

The  optimal  values  of  the  moael  constants  founa  in  this 
study  are  quite  different  from  the  values  previously  usea 
and  emphasize  the  need  for  careful  evaluation  of  the  con¬ 
stants  in  different  conditions.  The  values  of  these  con¬ 
stants  determined  by  Garwood  (1976)  were  ra-  =  7.6,  p_  *  m^/ra 
=  1,  r  =>0.5  and  O  »  0.1m  .  The  Garwood  model  has  previ¬ 
ously  been  tested  with  OS  Papa  data  by  Garwood  ano  Aaamec 
(1982)  ,  who  specified  m^  »  0.75,  p  =  1.0,  r  *  0.5,  ana  i^  - 
0.2m  .  The  ocean  thermal  structure  and  surface  flux  infor¬ 
mation  used  in  the  present  study  is  quite  aifferent  from 
that  used  in  the  OS  Papa  simulations.  The  ocean  thermal 
structure  used  in  the  OS  Papa  cases  was  aerivea  from  XBT  oo- 
servations  taken  at  the  same  site  approximately  every  two 
days.  The  thermal  structure  for  this  stuay  was  ootainea 
from  monthly  objective  analyses  of  XBT  ooservations  on  hori¬ 
zontal  surfaces  of  constant  depth.  Also,  the  XBT  observa¬ 
tions  can  resolve  temperature  gradients  with  Az  <  5m, 
whereas  the  TRANSPAC  analyses  ao  not  resolve  graaients  witn 
Az  <  20m.  In  the  OS  Papa  simulations  of  Garwooa  ana  Aaamec 
(1982) ,  the  surface  stress  and  the  heat  fluxes  were  calcu- 
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lated  from  3-h  observations  of  SST,  air  temperature,  oew 
point  temperature,  cloud  cover,  wind  speed  anci  aireccion. 
The  wind  speed  and  direction  were  aojusteo  to  10  m.  For 
this  study,  the  surface  stress  is  calculated  from  a  mooei 
produced  surface  wino  speeo  adjusteo  to  19  ra.  Tne  total 
beat  flux  and  solar  radiation  flux  are  calculateo  froai  ».iooei 
heat  and  moisture  fields. 

Model  sensitivity  to  variations  of  tne  values  of  tne 
model  constants  may  be  juuged  by  comparing  time  series  of 
SST  generated  using  different  values  of  the  constants  with  a 
time  series  of  the  analyzeo  SSTs.  The  improvement  in  mooel 
performance  using  the  values  oetermineo  in  tnis  stuay  com¬ 
pared  with  either  the  Garwooo  values  of  tne  constants  or  tne 
OS  Papa  values  is  illustrated  in  Figs.  3a,  o  ano  c.  In  Fig. 
3a,  using  values  which  are  typical  for  tnis  study  (m^  *  2.0, 
«  2.0,  r  ■  0.5,  if  ■  0.1  m”*) ,  the  mooel  accurately  simu¬ 
lates  the  annual  cycle  of  sea  surface  temperature.  Using 
the  OS  Papa  values  (m^-  0.75,  Pj*  1»0,  r*0.5,  if* 

0.2m**),  the  maximum  summer  temperature  is  coo  hign  aue  to 
too  little  wind  stirring  (m^  •  0.75)  and  aue  to  tne  snaliow- 
er  extinction  depth  {i(  ■  0.2ra*)  used  in  aosorpcion  of  tne 
solar  radiation  (Fig.  3b).  Further  mere  is  less  variaoili- 
ty  in  the  summertime  mixeo  layer  depth  ana  tne  rate  of  tne 
fall  mixed  layer  deepening  is  less  than  in  Fig.  3d.  In  Fig. 
3c,  using  the  values  aeterminea  oy  Garwooo  (m^  -  7.6,  p^  = 
1.0,  r  =  0.5,  if  *  0.1m**),  the  mooel  ooes  not  warm  enougn  in 
the  summer.  This  is  due  to  tne  fact  that  mj  is  too  large, 
causing  the  effect  of  the  wind  stirring  ro  oe  too  great  ana 
keeping  the  mixed  layer  mixed  too  oeep.  The  variaoiiicy  in 
the  summertime  mixed  layer  oepth  anu  tne  race  of  cne  rail 
deepening  is  much  greater  than  in  Fig.  3a.  Also  tne  winter 
mixed  layer  is  too  deep  compareo  to  Fig.  3a.  Anocner  impor¬ 
tant  difference  among  the  three  simulations  is  oeptn  of  tre 
summer  mixed  layer  and  the  time  of  the  occurence  of  cne 
spring  transition.  The  spring  transicion,  wnicn  is  cne 
abrupt  and  severe  retreat  of  the  mixeo  layer,  occurs  near 
the  end  of  February  for  the  simulations  using  tne  constancs 
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determined  in  this  study  and  the  OS  Papa  constants. 
ever,  in  the  simulation  using  the  Garwood  constancs,  tue 
transition  is  delayed  until  mid-April  and  the  upper  ocean 
never  achieves  the  shallow  mixed  layer  and  strongly  straci- 
fied  seasonal  thermocline  that  are  characteristic  or  tne 
summer  mixed  layer  in  mid-latitudes. 

5.  SUMMARY  AMD  CONCLUSIONS 

The  optimal  values  of  the  constants  for  any  turoulence 
model  are  expected  to  vary  for  different  specif icacions  of 
the  atmospheric  forcing  and  the  ocean  thermal  scructure.  iVe 
have  shown  that  the  values  of  the  mooel  constants  estimated 
from  laboratory  data«  and  Ocean  Weather  Snip  ooservacions 
yield  quite  different  results,  when  applied  to  the  same  sim¬ 
ulation  (Figs.  3a,  3b  and  3c) ,  than  the  "optimal”  values  de¬ 
termined  here.  Whether  it  would  be  better  to  alter  cne  val¬ 
ues  of  the  exchange  coefficients  for  the  surface  rluxes 
rather  than  the  values  of  the  closure  constants  is  largely  a 
philosophical  question.  Knowledge  of  surface  fluxes  ano 
geophysical  turbulence  is  not  yet  adequate  to  distinguish 
between  the  two  approaches.  We  nave  cnosen  the  second  ap¬ 
proach  largely  for  its  inherent  simplicity  ano  to  demon¬ 
strate  the  importance  of  understanding  the  potential  inter¬ 
actions  between  the  model  and  the  data  used  to  inicialize 
and  force  the  model. 

The  optimal  values  may  also  change  for  different  atmos¬ 
pheric  and  oceanic  regimes,  because  all  possible  pnysical 
processes  are  not  parameterized  in  the  model.  We  have  snown 
that  the  Garwood  model  is  relatively  insensitive  to  varia¬ 
tions  in  the  model  constants  within  a  fairly  broad  region  of 
parameter  space.  However,  the  year-to-year  variations  in 
the  optimal  values  of  the  model  constants  are  significant. 
In  particular,  1976  was  found  to  be  significantly  difrerenc 
than  1977  and  1978  by  all  three  statistical  tescs  employed 
in  this  study.  This  reflects  the  fact  that  the  atmospheric 
forcing  and  the  response  of  the  ocean  thermal  scructure  were 
quite  different  among  the  three  years  (White,  et.  al. ,  ?980; 
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Haney /  1980;  ana  Elsbercy^  1983)  .  An  iiuprovemenc  in  muoej. 
skill  could  be  achieved  by  calculating  the  optimal  values  oi 
the  model  constants  for  each  year.  Because  the  atmospneric 
forcing  and  the  ocean  thermal  scructure  are  calcuiateu  in 
the  same  manner  for  all  three  years,  the  variauions  in  cne 
optimal  values  of  the  constant.'^  probaoly  represent  tne  ef¬ 
fects  of  physics  not  incorporated  in  the  moael  ana  noise  in 
the  analyses  of  ocean  thermal  structure  ana  tne  surface 
fluxes.  For  these  reasons,  the  same  values  will  oe  usea  for 
all  experiments  using  this  type  of  forcing  in  the  North  Pa¬ 
cific.  The  values,  p  =  4.0,  ro_  =  2.0,  r  ®  0.5,  ana  iT  »  0.1 
•I  s  » 

m  ,  were  chosen  as  being  representative  or  the  swarms  of 
points  from  the  three  years  that  do  not  airier  significantly 
at  the  90%  confidence  level. 

Three  statistical  tools,  the  pool  ana  permutation  pro¬ 
cedure,  the  bootstrapping  technique  and  the  stuaent's  t- 
test,  were  used  to  determine  whether  aifferent  parameter 
values  gave  statistically  different  moael  results.  In  gen¬ 
eral,  the  bootstrapping  technique  generated  slightly  oroaaer 
reference  distributions;  however,  it  influencea  the  infer¬ 
ences  arawn  from  the  tests  in  only  one  case  (1976:  m,  *  3.5 
vs.  m^  s»  2.0).  The  student's  t-test,  on  the  otner  aana, 
lead  to  rejection  of  the  null  hypothesis  in  far  too  many 
cases.  This  points  out  the  nazard  of  assuming  Gaussian 
pdf's  for  this  type  of  data.  We  conclude  tnat  statistical 
tests  which  put  less  weight  on  the  inaepenuence  of  the  uata 
points  are  more  useful  for  stuaies  of  this  nature. 

The  absorption  of  solar  raaiation  with  depth  is  not 
considerea  one  of  the  funaamental  aspects  of  geopnysical 
turbulence  modeling.  However,  we  have  shown  tnat  its  proper 
specification  is  important  for  simulations  of  the  annual  cy¬ 
cle  and  crucial  to  simulations  of  the  summer  mixea  layer. 
The  parameterization  employed  in  the  Garwooa  moael  reproduc¬ 
es  the  results  of  more  complicateo  formulations  ana,  witn 
optimized  coefficients,  improves  the  summer  ana  annual  simu¬ 
lation  of  the  upper  ocean  temperature. 
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